Effective permittivity of mixtures of anisotropic particles 
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We use a new approach to derive dielectric mixing rules for macroscopically homogeneous and 
isotropic multicomponent mixtures of anisotropic inhomogeneous dielectric particles. Two factors 
of anisotropy are taken into account, the shape of the particles and anisotropy of the dielectric 
parameters of the particles' substances. Our approach is based upon the notion of macroscopic 
compact groups of particles and the procedure of averaging of the fields over volumes much greater 
than the typical scales of these groups. It enables us to effectively sum up the contributions from 
multiple interparticle reemission and short-range correlation effects, represented by all terms in the 
infinite iterative series for the electric field strength and induction. The expression for the effective 
permittivity can be given the form of the Lorentz-Lorenz type, which allows us to determine the 
effective polarizabilities of the particles in the mixture. These polarizabilities are found as integrals 
over the regions occupied by the particles and taken of explicit functions of the principal components 
of the permittivity tensors of the particles' substances and the permittivity of the host medium. 
The case of a mixture of particles of the ellipsoidal shape is considered in detail to exemplify the 
use of general formulas. As another example, Bruggeman-type formulas are derived under pertinent 
model assumptions. The ranges of validity of the results obtained are discussed as well. 

PACS numbers: 42.25.Dd, 77.22. Ch, 82.70.Dd, 82.70.Kj 
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I. INTRODUCTION 

The study of effective permittivity of heterogeneous 
systems holds an important position in various fields of 
physics and technics, for its results find wide-spread ap- 
plications in composite material engineering, biochemi- 
cal technology, medical diagnostics, etc. In theoretical 
research, the simplest and much discussed model con- 
siders a heterogeneous system as a mixture of fine par- 
ticles of the disperse phase embedded into a continu- 
ous host medium. The development of it began with 
the case of dilute mixtures of small spherical inclusions 
more than century agoi and has resulted, in particular, 
in the classical Maxwell-Garnett (MG) mixing rule and 
its various modifications (see Ref. for an analysis of 
relevant physical concepts and ideas from a historical 
perspective, and Refs. 3, 4 El for a review of major re- 
sults). It has been shown so far that: (1) The MG for- 
mula can incorporate multiple scattering effects (see, for 
instance, analytical results^ obtained in the quasistatic 
limit within a mean- field approximation). (2) For certain 
configurations of the disperse particles, it remains accu- 
rate for high-concentrated mixtures; ^ 10 ' 11 ! 12 in which 
strong electromagnetic interaction is significant and for 
which another - the Bruggeman mixing rulei^ - is often 
believed to be superior to the former. 

The case of mixtures of anisotropic particles remains 
little-investigated (see reviews^! after whose appearance 
the state of the art has not changed much). As far as 
we know, the existing attempts at taking the particles' 
anisotropy into account usually reduce to or heavily rely 
on different kinds of one-particle approximations, includ- 
ing their combinations (see, for instance, Refs. [T^fTol and 
references therein). A typical example of such approxi- 



mations is the use of the one-particle polarizability, de- 
scribing the response of a solitary particle to a uniform 
electric field, instead of the effective polarizability of the 
particle in the mixture. It is evident that such an ap- 
proach is tolerable only for diluted gases of anisotropic 
particles. In sufficiently concentrated mixtures, both 
multiple polarization effects and many-particle correla- 
tions in positions and orientations of the particles come 
into play. As a result, finding the effective polarizability 
becomes a many-particle problem, which is equivalent to 
the original problem of finding the effective permittiv- 
ity of the mixture. Correspondingly, neither the effective 
polarizability as a function of the dielectric and geomet- 
ric parameters of disperse particle, nor the interrelation 
of the two factors of anisotropy of the particles - non- 
sphericity of the shape and anisotropy of the substance 
- can be determined consistently within a one-particle 
approximation. Yet we are unaware of any practically 
important attempts at approaching these problems using 
the methods of multiple-scattering theory (in contrast to 
the case of spherical inclusions, for which see review^ 
key works^JI^^ and Refs. QjIJ. This fact is read- 
ily explained by the lack of knowledge of an infinite set 
of the correlation functions for concentrated systems of 
anisotropic particle. Even if such information were avail- 
able, practical calculations would be extremely difficult 
and would probably be limited to estimations of several 
corrections to the Born approximation. 

Recently^ we proposed a new approach to analysis of 
the long-wavelength value of the effective permittivity of 
finely dispersed mixtures. The idea was to avoid exces- 
sive theoretical refining on polarization and correlation 
processes that occur within the system on particle-size 
and interparticle-distance scales by averaging their con- 
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tributions out over macroscopic regions reproducing the 
properties of the entire system. The appropriate pro- 
cedure is based upon the notion of macroscopic com- 
pact groups of particles and the averaging^ 2 - of fields over 
volumes much greater than the typical scales of these 
groups. By applying it, we carried out^i a rigorous 
analysis of the effective permittivity of a concentrated 
mixture of spherically symmetric dielectric balls with 
piecewise-continuous radial permittivity profile. Later, 23 
the method was applied to systems comprising nonspher- 
ical inclusions with scalar permittivity. It was also shown 
that both the MG and Bruggeman mixing rules can be 
reconstructed with it. 

In the present report, approach 2 ^ is developed for 
macroscopically homogeneous and isotopic mixtures of 
anisotropic dielectric particles whose dielectric proper- 
ties are described by permittivity tensors and which are 
embedded in a host medium with constant scalar per- 
mittivity; the particles are assumed to be hard, measur- 
able, and, in general, inhomogeneous. It is shown that 
the averaged contributions from all-order reemission and 
short-range correlation effects within such a system can 
be effectively summed up. As a result, the effective static 
permittivity of the system is obtained as an explicit func- 
tion of the parameters of the model. It can be given the 
form of the Lorentz-Lorenz type, the polarization prop- 
erties of the particles being characterized by their effec- 
tive polarizabilities in the mixture. The latter are found 
via the geometric and dielectric parameters of the parti- 
cles and the host medium; as an example, a mixture of 
particles of the ellipsoidal shape is considered in detail. 
Finally, Bruggeman-type formulas are shown to follow 
from the general formulas under special choices of the ef- 
fective medium, and the ranges of validity of the results 
obtained are discussed. 



II. BASIC RELATIONS FOR ELECTRIC FIELD 
AND INDUCTION 

To begin with, we consider the problem on propagation 
of electromagnetic waves in a finely dispersed mixture 
with local permittivity ffjfe(r) = e^Sik + fe,fc(r). Here, Eq 
is the permittivity of the host medium and 8ei k (r) is the 
contribution caused by compact groups of the disperse 
anisotropic particles. By a compact group we understand 
any macroscopic region within which all interparticle dis- 
tances \r\ — r p | are small as compared to the wavelength 
of the probing wave in the mixture: y/eoko\r\ — r M | << 1, 
where fco is the wave vector of the wave in vacuum (see 
Fig. []}. According to Ref. Hi], we expect that it is multi- 
ple reemission and short-range correlation effects within 
such groups that form the effective permittivity of the 
mixture in the static limit y^fco ~~ * 0, where even quite 
large groups of disperse particles become compact. With 
respect to an external static field, such groups can actu- 
ally be treated as point-like; fluctuations of macroscop- 
ically large numbers of particles contained in them, as 
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FIG. 1: Homogenization of a mixture in terms of compact 
groups, viewed as macroscopic many-particle regions with lin- 
ear sizes d much less than the wavelength A of probing radi- 
ation. 



well as correlations between the groups, can be ignored. 
Correspondingly, the deviation Sei k (r) of the local per- 
mittivity in the mixture from the permittivity of the host 
due to the presence at the point r of a compact group of 
N identical disperse particles, occupying measurable re- 
gions fl\ with individual volumes v, can be modelled as 

N 

Se lk (r) = ^Ae lfe (r,fl A )n(r,fl A ), (1) 

A=l 

where A£^(r, Q\) — ei k (r,Q\) ~ eo^ife is the deviation 
of the local permittivity in the mixture from the permit- 
tivity of the host due to the presence of the Ath parti- 
cle at r, £jfc(r, stands for the permittivity tensor of 
the substance of this particle, n(r, is the character- 
istic function of the region fl\ [Il(r, = 1 if r G fl\ 
and LT(r, — otherwise], and 5i k is the Kronecker 
delta. In what follows, the tensors 5ei k (r), Si k (r, 
and Aeife(r, fl\) will also be designated as Se, s, and Ae 
respectively. 

The equation for the i component of the electric field 
of a wave in the mixture can be written as 2 ^ 

AEi + kfeoE, - V. t V fe S fc - -k 2 5e lk E k ; (2) 

henceforth, the dummy suffix convention applies. Eq. ([2]) 
is equivalent to the integral equation 

E t (r) = £7<°>(r) - k 2 ,J dr'T t] (r,r')8e jk (r')E k (r'), (3) 

v 

where (r) = Ef exp (i^/Foko ■ r) is the i component 
of the incident wave field in the mixture, E® is the 
i component of the field amplitude Eo, and Ty(r) = 
— (k^Sij + ViVj) exp (ik yfsQr) / AirkQEor are the com- 
ponents of the electromagnetic field propagator [Green's 
tensor for Eq. ©]; the integral in Eq. © is taken over 
the volume of the mixture. 
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With the iterative procedure applied, the solution of Our proof uses the representation 24 
Eq. can be represented in the form 



E i (r) = E^ (r)+^^(r) 



(4) 



E\ n \v) = {-kl) n J dv 1 J dv 2 ... J ' dr n T ij (T,T 1 )Se jk {T 1 ) 



V V 



xT k i(r 1 ,r 2 )Sei m (r 2 ) . . . T rs (r n -i, r„)& st (r„) J B t (0) (r n ). 



The i component of the electric induction vector at the 
point r in the mixture is given by 



Di(r) = [saSik + 6e ik (r)} E k (v). 



(6) 



Assuming that the mixture as a whole is macroscop- 
ically homogenous and isotropic, we define its effective 
permittivity e e ff in the standard way— as the propor- 
tionality coefficient in the relation 



Di — e c ff E. 



(7) 



where the bars indicate averaging by the rule Ai — 
y J Ai(r)dr; the volume V is implied to be much greater 

than the volumes of compact groups. Following Ref. l22j, 
we accept that for finely dispersed mixtures, the values Ai 
are equal to the corresponding statistical averages (Aj(r)) 
over the positions and orientations of disperse particle. 

Now, under the suggestions made, we can prove quan- 
titatively the following two statements: (f) In the lim- 
iting case y^fco — * 0, the values Ei and Di are deter- 
mined only by the multiple reemission and short-range 
correlation effects inside compact groups. (2) The cor- 
responding all-order contributions to e c g can be singled 
out from the iterative series by formal replacing the fac- 
tors kyTij(r) in the integrals for Ei and Di with the 
expressions 5(r)Sij /3sq, where <5(r) is the Dirac delta 
function. This replacement simply reflects the fact that 
within a macroscopic approach, the specified contribu- 
tions are formed by those ranges of the integration vari- 
ables where the electromagnetic field propagators reveal 
a singular behavior. 



3^M(r) + 2 ^(^-3e^) (8) 



for the propagator 2y(r), which is valid (see Ref. [25| for 
mathematical details) on a set of scalar, compactly sup- 
ported, and bounded functions #£ Q| a(r) in the sense that 
(5) there holds the relation 



lim / drfcgTy(r)<fe a/3 (r) = lim / drfc 2 ^ (r)fc Q/3 (r). 

V V 

In Eq. (|8|) , e$ stands for the i component of the unit vector 
e = r/r; for brevity, the first and the second terms to the 
right of the equality sign will also be denoted by f W and 



-(2) 



respectively. 



With the aid of ((8]), the contribution to the statisti- 
cal average (Ei(r)} from the nth iterative step can be 
represented in the static limit y^fco ~~ * as 



(n) 



E 



The addend 



(r)) = (^(r) 



E 



(n;l,2) 



(r) 



(9) 



E 



(«;!)/ 



(r)) = y-^-J (5e ik (r)Se km (v) . ..Se rt (r))E° 

(10) 

represents the statistical average of a product of n fac- 
tors 5e related to a single compact group; it is obtained 
by replacing all factors k^T in (0 by the their most sin- 
gular parts f^ 1 ). The addend (E^ n '' 1 ' 2 \r)^ is the sum 

of all integrals containing at least one factor under 
the integral sign. Each expression containing m < n fac- 
tors of the f ( 2 ) type results from taking (n — m) integrals 
with factors in their integrands; it is therefore pro- 
portional to the expression 



/ drj dr 2 ...f dr m (r - n) f ^ (n - r 2 ) . . . (r ro _i - r m ) 

(fc(n)) n (fc(r 2 )) 12 . . . (fe(r m )) 4m )E°, 



V V V 

X 



where ii+i 2 + ... + i m = n. Since the spacial and orienta- tional correlations between macroscopic compact groups 
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in a macroscopically homogeneous and isotropic mixture 
are negligibly small, the many-point (involving different 
compact groups) correlator in this expression can be fac- 
torized as a product of one-point correlators, each of 
which is related to a single compact group. The one- 
point correlators are equal to linear combinations of ex- 
pressions constructed of Kronecker deltas and satisfying 
certain symmetry constraints on permutations of their 
indices. In particular, the average (5sj k (ri)) = ASj k is 
symmetric with respect to permutations j <-> k, the aver- 
age (fe jfe (r 1 )fc im (r 1 )) = B6j k 6 lm + C (SjtSkm, + 6 jm 5 kl ) 
is symmetric with respect to permutations j «-> k, I «-+ m 
of individual indices and permutations (J, k) «-> (l,m) of 
their pairs, etc. The coefficients A, B, C, . . . in these rela- 
tions depend on the physical parameters of the host and 
disperse particles, but not on the coordinates of the rel- 
evant compact group. As a consequence, the integrals in 



Note that there are, respectively, a and a + 1 factors 
under the bars in the sum from 1 to oo and that from 
to oo. 



III. EFFECTIVE PERMITTIVITY OF 
MATRIX-PARTICLE MIXTURES 

For a macroscopically homogeneous and isotropic mix- 
ture, the averages Ei and Di are proportional to the cor- 
responding component Ef of the external field. Since 
different terms in series (|13[) and (|14|) are independent, it 
is reasonable to suggest that they have the structure 



feife(r)fe fem (r) . . . fe r t(r) = A a 8 lU (15) 

where the subscript a specifies the number of factors un- 
der the bar. The summation over the indices i — t yields 



A„ = -Tr (fc(r))". (16) 



^-Ej n ' 1,2 ^(r)\ reduce to those taken of factors alone. 

Taking into account the explicit form of p 3 ', we see that 
these integrals vanish after integration with respect to 
the angles. Thus, 

(if ;1 ' 2) (r)) = 0. (11) 

In view of the equality of the averages (-Ei(r)) and Ei, 
Eqs. ©-d) give 

Et\r)= (-3^) 5e lk {v)8e km {r)...5e rt {r)E t . (12) 

Similar reasoning is used to analyze the averages 
(Dj(r)) and Di. From Eqs. dH)-©, we finally obtain: 



(13) 



(14) 



This trace is easily found, for in the typical expression 



5e lk (r)Se km (r) . ..5e rt (r) 
n n N 1 r 

= E E ••• E vj dr ^ik(r,n Xl )n(T,n Xl ) 

Al — 1 A2 — 1 A CT — 1 y 

x Ae km (r , n X2 )II(r, Q^)... Ae rt (r, n Xtr )H(r, n K ) 

all addends with two and more differing values of the 
indices Ai, A2, . . . , A CT are zero (the regions occupied by 
hard particles never overlap). If all of these values are 
equal, then the relation [II (r, ^a)] ct = n(r, f2^) (er any 
natural number) gives 



Tr (fc(r)) CT = E^ / drAe lk (r,n x )Ae km (v,n x )... 

xAs ri (r,Q x )=nJ drTr (Ae(r, . 

V 

Here, n — N/V is the particle concentration, the last 
integral is taken over the region occupied by a single 
particle, and the integrands are assumed to be wise- 



S it 



E 

cr=l 



3e 



S e lk {r)5e km (v) . . . 5e rt {v) 



D, 



00 / 1 \ CT 

£oSu + £oE ( "3^ ) Se lk (r)Se km (r) . . . 6e rt (r) 



00 / 1 \ a 

+ E ( _ 3^) fc ife( r ) fc fcm(r) • ■ ■ Se pr (r)ds r t(r) 



E-. 
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continuous. Thus, 

A<, = ~J rfrTr (Ar(r.<>\ 



or, in terms of the principal components Aeu(r) of Ae 
(t = l,2, 3), 

A ff = | / dr{[A £ll (r)] CT + [Ae 22 (r)] CT + [Ae 33 (r)y }, a > 1. 

(17) 

It follows from Eq. (jTTJ) that the field |13]) and in- 
duction (fT4|) are represented by infinite geometric series. 
Summing them up and using definition (JT]), we obtain a 
formula of the Lorentz-Lorenz type: 



, , 8tt \ / 4tt 
£ c ff = £o I 1 + -^-^a c ff 111 — — na e g 



where 



"eff = - (an + a-22 + "33) j 



An J 2s 



( r ) ~ £0 



i = l,2, 3. 



(18) 



(19) 



(20) 



According to Eqs. (|18p - (l2"01) , the effective polarization 
properties of an anisotropic particle in a finely dispersed 
mixture are described by the quantity a c s. The latter 
can be treated as the effective polarizability of the par- 
ticle in the mixture. The value of a e ff is found as the 
arithmetic mean (| 19|) of the quantities an , related to the 
corresponding principal components En (r) of the permit- 
tivity tensor of the particle's substance. However, it is 
physically incorrect to interpret an as the principal com- 
ponents of a certain tensor, which could be called the 
effective polarizability tensor of disperse particles. The 
effective polarizability a c s is in fact a scalar quantity, 
contributed to by different physical mechanisms; their 
combined effect is given by (TT9|) . 

Generalization of Eqs. (|18p - ([20|) to the case where 
a mixture comprises particles of different sorts a = 
1, 2, . . . , s, with individual volumes v a and concentrations 
n a , is evident: 



£cff = £0 ( 1 + Uaa °s J ( 1 ~ naQ £ 



eff i 



a { a 1 a , a \ 

"off = 3 ("11 + "22 + «33J > 



(21) 
(22) 



AttJ 2e + t 



£0 



?i(r)' 



1, 2, 3, (23) 



where e H (r) are the principal components of the permit- 
tivity tensor of the substance of particles of sort a, and 
a" 9 can be interpreted as the effective polarizabilities of 
these particles. 

For particles filled with homogeneous anisotropic di- 
electrics, Eqs. (|20|) and ([23]) take the form (no summation 
over a) 



_2_ v £ u - £ o 
4tt a 2e + ef, 



(24) 



It is interesting to note that quantities (j24|) for an 
anisotropic particle are formally equal to the principal 
components of the one-particle polarizability tensor for a 
ball filled with the same dielectric and having the same 
individual volume. 

The effective permittivity of a mixture of homogeneous 
anisotropic particles is 



£ c ff 



£0 



3^'"^ 2e 

a=l i=\ u 



/ I" 3 £° - £ \ 1 



(25) 



where c a is the volume concentration (fraction) of parti- 
cles of sort a. 

For two-component mixtures (s = 1), the result (|25j) 
agrees with some of the rules known in the literature. 
Two particular examples are of interest. (1) The disperse 
particles consist of a substance with isotropic dielectric 
properties (en = £22 = £33 = £)■ Then Eq. ([23]) re- 
duces to the classical MG mixing rule, no matter what 
the shapes of the particles are: 



£cff = £0 + 3C£ 



£0 



+ 2e - c(e - e ) ' 



(26) 



(2) The disperse particles consist of a uniaxial sub- 
stance, with one permittivity value (£33 = £11) in one 
preferred direction and another in all perpendicular di- 
rections (£n = £22 = £_l). Now, Eq. (f25|) takes the form 

£ ff = £0 + 3c£ 

(£_L + 2£p)(£|| - £p) - 2£p(£|| - £j_) 

X (1 - c)(£L + 2£ )(£|| + 2£ ) + C£ (£ ± + 2ejj + 6£ ) ' 

(27) 

Eq. (|27p first appeared in Ref. [IJ, where an anisotropic 
version of the MG approximation was developed for the 
special case of a mixture of randomly-oriented uniaxial 
spherical particles. 

It should be emphasized that the mixing rules (|18[) . 
([2"Tjl . and (|2"5)) - (|2"T|) were obtained for sufficiently con- 
centrated mixtures of anisotropic particles. In deriving 
them, no constraints on the value of the difference be- 
tween the permittivities of the particles and host were im- 
posed. To evaluate feasible restrictions on concentration 
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values, we note that these rules were derived for macro- 
scopically homogeneous and isotropic mixtures with the 
local permittivity tensor given by Eq. ([1]). The model 
(fTj) assumes that the particles of a mixture retain their 
entities to the greatest extent possible; in other words, 
they are viewed as clearly distinct inclusions embedded 
into a host medium (matrix). For spherical particles 
and under certain restrictions (see recent experimental, 11 
numerical^ and theoretical^ estimates), the classical 
MG mixing rule can be accurate for concentration val- 
ues as high as c^x ~ (0-4 ^ 0.5). In the case of non- 
spherical particles, the shape effect comes into play: for 
concentrations exceeding a certain value c max , different 
orientations of elongated particles become hindered by 
the neighboring particles; as a result, the individualities 
of such particles start being disturbed. Given the great- 
est linear size, G? max , of a particle, c max can be estimated 
as c max sa cfP ax i>/u maX5 where v is the volume of the par- 
ticle and w max = nd^^/f) is the minimum volume of 
an imaginary sphere admitting of free rotations of this 
particle. Indeed, if the bulk of each such sphere were 
completely filled with the substance of the particles, the 
above-mentioned results for mixtures of spherical parti- 
cle would apply. However, the actual concentration of 
the filling substance is w max /u times less than cfP^ x . 

We suggest that the upper bound of the validity range 
in concentration for Eqs. ([18]). (|2Tj) . and d25j)— (J2TJ) is no 
less than c max . For certain packings of particles, even 
relatively long, and also for mixtures with wide spreads 
in the shapes and sizes of particles, the value of the this 
bound is expected to increase. 



A. Mixtures of ellipsoidal particles 

Consider in detail a mixture of ellipsoids with semiaxes 
di, d 2 , d 3 ; the ellipsoids are filled with a homogeneous 
anisotropic dielectric with permittivity tensor Sik- Be- 
sides practical applications, this case attracts interest as 
an example where the effective polarizability a e g can be 
expressed explicitly through the geometrical parameters 
and the principle components of the one-particle polar- 
izability tensor of a solitary ellipsoid. These components 
can be calculated theoretically and measured experimen- 
tally. 

For the sake of simplicity, let us suppose that the prin- 
cipal axes of the tensor Sik coincide with the principal 
axes of all ellipsoids. The principal components of the 
one-particle polarizability tensor of such an ellipsoid in a 
host medium of permittivity Sq are^ 2 - 



1 

7« = 7Z V 



£o 



47r s + Li(su - s ) 



i = 1,2,3, 



(28) 



where Sa are the principal components of s^ and are 
the depolarization factors. The latter are given by the 



integrals 



did 2 d 3 



dit 



2 J {u + dl)R{u) 
o 



R(u) = ^{u + dl){u + dl){ U + df), 



and satisfy the relation 



Li + L 2 + L 3 = 1. 



(29) 



Using Eqs. ([24]), ([28]), and ([29]), the desired relation is 
easy to find: 



1 3 

« off = 



l + (l-3ii) 



7i-. 



did 2 d 3 



(30) 



If there are different sorts of ellipsoidal particles in a mix- 
ture, Eq. ([3TJ]) and its particular versions (pTT]) . ([3"3"]) hold 
for each sort separately. 

Eq. ([30]) reveals that the effects of two factors of 
anisotropy, the shape of the particles and anisotropy of 
the dielectric parameters of the particles' substance, on 
a e g are interlinked in an intricate way. In the case of 
a mixture of anisotropic balls, which is free of the shape 
anisotropy effects and where Li — 1/3, Eq. ([30]) takes the 
form 



1 



(ill + 722 + 733) ■ 



(31) 



If the ellipsoids are filled with a homogeneous isotropic 
dielectric with scalar permittivity e, then the principal 
components of their one-particle polarizability tensor are 



0ii = ~r v 



£() 



i = l,2,3, 



47r " so + Li(s - so) 
and their effective polarizability in a mixture is2£ 



,1 1 1 
c*cff = 3 ( — + — + — 
,Pn P22 P33 



(32) 



(33) 



It follows that the shape anisotropy results in a nonlinear 
relationship between the effective polarizability and the 
principal components of the one-particle polarizability. 
For a mixture of homogeneous isotropic balls, the values 
of these polarizabilities are equal. 

It should be remembered that results ([30]) . ([3T]) , and 
(|3"3"]) are based on the mixing rule ([2l>]) . Except for the 
case of spherical particles, its functional form is signif- 
icantly different from those typical of one-particle con- 
siderations. For instance, one of the most popular MG- 
type formulas for a two-component mixture of randomly- 
oriented homogeneous ellipsoids reads^ 



1 3 

£eff = £0 + Z c£ 0^2 



£0 



1-ic- 



—j 1 £ + Li(s - s ) 
Li(s - s ) 



3 ^ so + L. L (s 

I— 1 



(34) 
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FIG. 2: S as a function of c for mixtures of (a) prolate 
spheroids (aspect ratio 2) and (b) oblate spheroids (aspect 
ratio 0.5) embedded in a uniform isotropic host. The permit- 
tivities eo and e are: 2 and 10 (curve 1); 2 and 100 (2); 2 and 
1000 (3); 10 and 2 (4); 1000 and 2 (5). 



A question arises how close are the predictions made 
by Eq. (|26|) [the particular case of Eq. (|25|) for two- 
component mixtures] and by other mixing rules, such as 
Eq. 

To evaluate the numerical discrepancies between 
Eqs. (|26[) and l|34|) , we have analyzed the relative dif- 
ference 6 = [(ejjp - eJP)/eJP] x 100% in the val- 
ues that they give for e c g. Fig. [5] illustrates the be- 
havior of 6 for mixtures of prolate {d\ — g?2 < d%) 
and oblate (d\ = c?2 > da) spheroids with aspect ra- 
tios d 3 /di = 2 (c max ps 0.125) and 0.5 (c max pa 0.25), 
respectively. Several facts are clearly seen for these sys- 
tems. (1) In the limiting case of diluted mixtures, S — ► 0. 
(2) Even if the contrast k = e/eq is extremely high, \6\ 
does not exceed 3 -j- 4% for c < 0.5c max and 5 -j- 7% for 
c £ (0.5c max , c max ). (3) For moderate values of c, |5| in- 
creases when either c, or k, or they both increase; depend- 
ing on k, \S\ can be as high as 10-^15%. (4) Compared to 
Eq. (f3"4"| . Eq. (f2T))) predicts a smaller magnitude for the 
ellipsoidal- sha pe influenc e on the formation of e c s- fo r a 

given c, — £o| < kffP^ — £ o|- The physical expla- 

nation for these tendencies is, in our opinion, as follows: 
multiple mutual polarizations and short-range correla- 
tions between ellipsoidal particles effectively reduce the 
anisotropy effects induced by a single particle. Unfortu- 
nately, the lack of numerical data for the systems under 
consideration does not allow us to check directly which 
of Eqs. (USD and jSJ) does better. 



IV. BRUGGEMAN-TYPE MIXING FORMULAS 



As the concentration of anisotropic particles in a mix- 
ture becomes greater than c max , two kinds of phenomena 
can be expected. On one hand, strong mutual influences 
can trigger various physicochemical processes, affecting 
the textures and dielectric properties of the particles and 
host medium. On the other hand, interparticle contacts 
can cause the formation of coarse aggregates of particles; 
as a result, asymmetry between the disperse phase and 
the host medium fades away. In the limit, the original 
mixture, including particles of s sorts, can be viewed as 
a system of s + 1 components, filling chaotically located 
regions with irregular complex boundaries. 

It was already mentioned that for concentrated mix- 
tures, the Bruggeman mixing rule is often believed to be 
superior to the MG-type mixing rules. However, the ac- 
curacy and range of validity of this rule have remained 
disputable (for relevant details, see extensive review^). 
To a great extent, the situation is explained by the lack 
of accurate electromagnetic solutions for disordered sys- 
tems; as a result, most arguments for the Bruggeman 
rule are based on a dipolar approximation, valid for low 
concentrations of particles. The purpose of this section 
is to shed some light on theses issues. Namely, we show 
that the Bruggeman-type formulas are obtainable within 
the proposed compact-group approach under appropri- 
ate assumptions. In addition, system-dependent devi- 
ations from the classical mixing rules are expected to 
occur, provided physical and chemical processes change 
the dielectric and structural parameters of the mixture 
constituents. 

To begin with, we remind that Eq. (fl"3"|) and (Q3J) are 
valid for any mixture satisfying the conditions of macro- 
scopical homogeneity and isotropy. They form the basis 
for analysis of particular structural models of mixtures, 
which arc formulated in terms of the local permittivity 
tensor £j/s(r). For instance, the above MG-type mixing 
rules (see Section IIII[) are based on the model expres- 
sion JTJ). Now, let us suppose that all components of 
a mixture, including the host medium, can be treated 
symmetrically in the following sense: there exists a ficti- 
tious (effective) medium, of permittivity e, such that the 
real mixture is macroscopically equivalent, in its dielec- 
tric properties, to a macroscopically homogeneous and 
isotropic system of nonoverlapping regions occupied by 
the components of the real mixture and embedded into 
this fictitious medium (see Fig. [3]). If combining the 
components into a real mixture does not change their 
properties, then the deviations of the local permittiv- 
ity in the system from the value e can be written as 
Ae ik {r,tt\ a ) = £ ik (r,Q,\ a ) - eS lk for regions Q\ a , occu- 
pied by particles of sort a, and as (eo — e)Sik for the rest of 
space, occupied by the host. Correspondingly, the local 
permittivity deviations Seikij) from e due to a compact 
group, of the above regions, present at the point r in the 
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FIG. 3: Effective medium approximation for the mixture FIG. 4: Modelling the mixture shown in Fig. \T\ for the case of 
shown in Fig. \T\ high particle concentrations. 



system can be modelled as 

s N a 

Se ik (r) = (e -e)S lk U(r)+^2 E ^aJ H(r, fi A J, 

a=l A a = l 

(35) 

where 

s N a 

n(r) = i-£ E n ( r '°Aj 

a=l A a =l 

is the characteristic function for the space occupied by 
the host of the real mixture. Carrying out manipula- 
tions similar to those in Section IIIII and introducing the 
total volume concentration c = X)a=i c a °f ^ ne disperse 
particles, we find: 

2e + £ cff 2e + £o 3 ' 

a— 1 

= 5(011 + 2^2 + Sis) i (37) 



Under the additional assumption that the fictitious ef- 
fective medium is macroscopically equivalent, in its di- 
electric properties, to the mixture itself, e = s c g, the 
left-hand sides in Eqs. (|36|) . (j39|) . and ([40]) vanish, and 
we arrive at generalizations of the classical Bruggeman 
fixing rule. Note also that particular forms of Eq. (|40l) 
are known in the literature (see, for instance, Ref. |26| . 
where it was obtained within a dipolar analysis for the 
case of identical spherical inclusions). Finally, for e = Eq, 
the results of Section IIIII are obtained. 

For sufficiently concentrated mixtures, large aggre- 
gates of particles can themselves be treated as structural 
units, with scalar (on the average) permittivities. Mod- 
elling such a mixture as a system of nonoverlapping re- 
gions occupied by substances (the host medium of the 
real mixture among them) with different permittivities 
and embedded into a fictitious medium of permittivity 
e (see Fig. @|, we again arrive at Eqs. (|3"6"|) - ([3"5)) . with 
e»(r) = e a (r), and Eq. (gDJ), with e a (r) = e a = const. 
Now, £ a (r) stands for the averaged-permittivity distribu- 
tion within the region of sort a, n a is the concentration 
of such regions, v a indicates integration over the region 
of sort a, and c a stands for the volume concentration of 
the substance of sort a. 



5S = A/^(r) 



4nJ 2e + eS(r)' 



i=l,2, 3. 



(38) 



CONCLUSION 



The quantities cf" ff can be interpreted as the effective po- 
larizabilities of particles of sort a in the fictitious effective 
medium of permittivity e. For mixtures of particles filled 
with homogeneous anisotropic or homogeneous isotropic 
dielectrics, we respectively have: 



£efl — e 
2e + £ ff 



(1-c) 



e - e 



1 



E c -] 



E £ ii 



£cff — e /-, n £0 
= (1 - c) 



2e + £ c ff 



a— 1 i—l 



E Ca < 



2e + £ ^ 2e + £ a 

a— 1 



(40) 



In the above formulas, the permittivity e is unknown 
and should be considered as an adjustable parameter. 



The main results and conclusions of this paper are as 
follows. 

(1) The compact-group method for study of effective 
dielectric properties of finely dispersed mixtures is gen- 
eralized to the case where the mixtures are composed of 
anisotropic inhomogeneous particle. Based on this gen- 
eralization, new mixing rules for mixtures of hard dielec- 
tric particles are derived and their ranges of validity are 
discussed; in appropriate cases, the rules obtained agree 
with the classical Maxwell-Garnet and Bruggeman rules. 
It should be emphasized that within our method, the 
contributions from multiple polarization and short-range 
correlation effects of all orders are effectively taken into 
account. 

(2) That fact that mixing rules of both types are ob- 
tainable (under appropriate approximations) within a 
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single formalism is taken by us as strong evidence of 
the validity of basic relations fl~3")) and (|14[) . If this is 
really the case, much emphasis in further effective per- 
mittivity studies should be put on modelling, in terms 
of the permittivity distribution, of the structure and di- 
electric properties of the constituents of a real mixture. 
Deviations of these properties from those of isolated con- 
stituents should manifest themselves as deviations from 
the classical mixing rules and their modifications ob- 
tained in the present work. 

(3) The averages in Eqs. (fTB")) and ([TJJ are in fact sta- 
tistical, according to the analysis in Section [II] Con- 
sequently, Eqs. (fTo) and (fT4|) can be used to develop a 



statistical theory of effective permittivity. In particular, 
proceeding in this way and using the general properties 
of particle distribution functions for macroscopically ho- 
mogeneous and isotropic systems of balls, one can infer 
that it is the hardness of the balls that lies at the heart 
of the classical MG mixing rule. 
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